Time-Reversal Symmetry-Breaking Superconductivity in Heavy Fermion PrOs 4 Sbi 2 

detected by Muon Spin Relaxation 
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We report on muon spin relaxation measurements of the 4/ 2 -based heavy-fermion superconductor 
filled-skutterudite PrOs4Sbi2. The results reveal the spontaneous appearance of static internal 
magnetic fields below the superconducting transition temperature, providing unambiguous evidence 
for the breaking of time-reversal symmetry in the superconducting state. A discussion is made on 
which of the spin or orbital component of Cooper pairs carries a nonzero momentum. 

PACS numbers: 74.70.Tx, 76.75.+i, 74.70.Dd, 74.25.Ha 



Many unconventional superconducting (SC) states, as 
in Ce- and U-based heavy-fermion (HF) compounds or 
high-T c cuprates, appear in close proximity to magnetic 
instabilities when a certain parameter (pressure, atomic 
doping, or oxygen content) is varied 0,13, El This fact 
strongly suggests that the attractive interactions binding 
electrons into Cooper pairs are mediated by magnetic 
moment fluctuations. As another possible pairing glue, 
fluctuations of quadrupole moments — distorted shapes of 
the electronic clouds of ions, are theoretically considered 
to be possible 0, Q . An interesting question to be ad- 
dressed is what is the nature of superconductivity in such 
a system. For this study, Pr-based compounds with a 
4/ 2 configuration are likely candidates, since nonmag- 
netic but quadrupolar active low-energy levels can be re- 
alized due to the crystalline-electric-field (CEF) effect; in 
5/ systems, CEF levels are less clear due to the tendency 
to be itinerant. 

One promising candidate material for this study is the 
recently found superconductor PrOs4Sbi2 6], which is to 
date the only known Pr-based HF superconductor, with a 
superconducting transition temperature of 1.82 K (here- 
after referred to as T& ) ■ The estimated electronic specific 
heat coefficient 7 = 350 - 700 mJ/K 2 mol 0,0 and the 
enhanced cyclotron-effective masses |SJ reflect the exis- 
tence of strong electron correlations. Specific heat (C), 
magnetic susceptibility (x), and inelastic neutron stud- 
ies provide evidence that PrOs4Sbi2 has a nonmagnetic 
ground state and a magnetic triplet excited state sep- 
arated by AE'cefAb = 8 K [ZL 0], which is 5 times 
larger than T c ±. In the temperature- versus-magnetic- 
field (T-vs-H) phase diagram, a field-induced ordered 
phase (hqH >4T) 9| appears close to the superconduct- 
ing phase (the upper critical field /iqH C 2 = 2.2 T). It was 
recently proven to be an antiferro-quadrupolar ordered 
phase by elastic neutron scattering measurements 10]. 
This fact strongly indicates that quantum quadrupole 
fluctuations of the Pr ions play an important role in 
realizing the HF superconductivity in PrOs4Sbi2, con- 



sidering that the T-vs-H phase diagram is analogous to 
those for the HF and cuprate systems, where a magnet- 
ically ordered phase exists close to the SC phase in the 
T-vs-pressure, -atomic-doping, or -oxygen-content phase 
diagram. This scenario is further supported by the en- 
hanced T c \ compared to T r = 0.74 K for a 4/° reference 
compound LaOs4Sbi2 |a.Hl). 

The remarkable unconventional SC properties de- 
scribed below suggest non-s-wave pairing in PrOs4Sbi2- 
Thermal conductivity measurements in magnetic fields 
rotated relative to the crystal axes indicate the pres- 
ence of two distinct superconducting phases with differ- 
ent symmetries with point nodes [l2j. The specific heat 
also suggests an additional phase transition at T C 2 — 
1.65 - 1.75 K < T c i EM Sb-nuclear-quadrupole- 
resonance (NQR) spin-lattice- relaxation rate 1/Ti ex- 
hibits no Hebel-Slichter peak just below T c i 0|. In con- 
trast, exponential T dependences of 1/Ti in the lower T 
region [llj and of the magnetic penetration depth deter- 
mined by a transverse-field muon-spin-rotation study [l4j 
indicate, however, an isotropic energy gap. These behav- 
iors, making it difficult to construct a simple picture for 
the SC state, suggest that a novel type of HF supercon- 
ductivity is realized in this compound. 

One important aspect for characterizing a SC state 
is whether time-reversal symmetry (TRS) is broken or 
not ^5|. In a state with broken TRS, the magnetic mo- 
ments of Cooper pairs are nonzero and align locally, and 
thereby a spontaneous but extremely small internal mag- 
netic field can appear. To detect such fields, zero-field 
(ZF) muon spin relaxation (/iSR) is the most powerful 
method, because of its high sensitivity, as small as 10 
fiT (TallTj. TRS-broken SC states are quite rare, since 
an unambiguous observation of such fields has been made 
only in Sr 2 Ru04 0|. In the low-T phase of Th-doped 
UBei3, the possibility of a magnetic phase coexisting 
with the SC state cannot be ruled out [lj| and, on an ob- 
servation in UPt3 , doubt has been cast later 0, E^ ■ 

In this Letter, we report on ZF-/iSR measurements for 
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FIG. 1: Zero-field /iSR spectra in PrOs4Sbi2 measured across 
the superconducting transition temperature T c i = 1.82 K. 
The corrected asymmetry was obtained by subtracting the 
temperature-independent contribution from the silver plate 
(0.057). The curves are fits to the model I given by Eq. (1). 
The relaxation rate becomes stronger below T c \. 



PrOs4St>i 2 . The results unambiguously reveal the spon- 
taneous appearance of an internal magnetic field in the 
SC state, providing clear evidence for broken TRS. This 
is the first observation of such superconductivity in un- 
doped intcrmctallic compounds. 

Conventional /iSR measurements were carried out 
down to 20 mK using a top-loading dilution refrigerator 
at the 7rA-port of the Meson Science Laboratory, KEK- 
MSL, Japan. Small single crystals of PrOs4Sbi2 with 
sizes of < 1 mm grown by the Sb-flux method [8j using 
raw materials of 4N(99.99% pure)-Pr, 3N-Os, and 6N-Sb 
were used. Clear de Haas-van Alphen (dHvA) oscilla- 
tions observed in one of the crystals are indicative 
of their high quality. The single crystals with random 
orientation were glued onto a 4N-silver holder covering 
~ 25 mmc6. The silver holder was held in a vacuum 
space of the refrigerator by screwing it to a copper rod 
thermally connected to the mixing chamber. Stray fields 
at the sample position were cancelled to within 1 fiT in 
all directions by using three pairs of correction coils and 
a flux-gate magnetometer. We implanted spin-polarized 
positive muons into the sample and the time evo- 

lution of the /i + spin polarization was measured using 
positron detectors positioned parallel and antiparallel to 
the initial polarization direction. 

The time evolution of the ZF muon spin polarization is 
shown in Fig.[TJfor T = 4.2, 2.0, and 0.02 K, where a non- 
relaxing background signal originating from /i + stopping 
mainly in the silver holder had already been subtracted. 
No signature of precession is visible, ruling out the pres- 
ence of a sufficiently large internal magnetic field as seen 
in magnetically ordered compounds. It is obvious from 
Fig-fflthat the- relaxation becomes stronger with decreas- 
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FIG. 2: Temperature dependences of the width of the internal 
field distribution A, the relaxation rate A in zero-field (ZF) 
and in ^io-Hlf = 0.01 T, and the specific heat divided by 
temperature C/T. The vertical lines indicate two consecutive 
transitions at T c i and Tci determined from the C/T data. 



ing T. We tried several functional forms to reproduce the 
ZF spectra, and excellent fits to the data were obtained 
using two models, I and IT The fitting using model I, 

KT/ 



P lt = exp(-At)Gf 1 (A,t), 
Gf T (A,t) = 1 + 1(1- A 2 t 2 ) exp(-iA^), (1) 



which is commonly used [TJ, [13, |2(j , is shown by the solid 
curves in Fig. □ The Kubo-Toyabe (KT) function Gf T 
is attributed to the muon-spin relaxation due to static 
randomly-oriented local fields at the muon site (caused 
by nearby nuclei in the normal state), and (A/7^) 2 
(7^ = 8.516 x 10 8 s _1 T _1 : the /i + gyromagnetic ratio) 
represents the second moment of the field distribution. 
The necessity of exp(— At) indicates the presence of an 
additional relaxation process. The best-fit values of A 
and A are shown in Fig. [21 along with the specific heat 
divided by temperature C/T measured using a single 
crystal from the same batch. Around the SC transition, 
C/T shows two types of anomalies, i.e., a sharp jump at 
T c i = 1.82 K and a kink structure at T C 2 — 1.64 K, in a 
slightly different way from the reported data 0, 0] . 

It is remarkable that A shows a significant increase 
with an onset temperature of around T cl , indicating 
the appearance of a spontaneous internal field correlated 
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with the superconductivity; this behavior is in marked 
contrast with that of A, which increases gradually with 
decreasing T across T cl and T c2 without showing any 
distinct anomalies. This observation provides unam- 
biguous evidence that TRS is broken in the SC state 
of PrOs4Sbi2. Note that the observed increase in A is 
smaller than the experimental uncertainty in the previ- 
ous ZF-/kSR measurement 14] . Provided that the in- 
crease is of electronic origin, the electronic (A e ) and nu- 
clear (A n ) contributions to A must be uncorrelated and 
thereby their respective contributions add in quadrature: 
A e (T) 2 + A 2 n = A(T) 2 . FromFig.El A„/ 7m = 1.7 x 10" 4 
T and A e (T -> 0)/^ = 1.2 x 1CT 4 T are derived. The 
temperature dependence of A e below T c \ agrees with that 
of the BCS order parameter Abcs(T), as is evident in 
Fig. |21 supporting our interpretation. 
The equally probable model II, 

P M = exp(-At)[~ + ^cos(^)]Gf T (A„,i), (2) 

also fits the data, yielding almost the same \ 2 statis- 
tics as model I; i.e., the two models are experimentally 
undistinguishable. This fact points to an alternative pos- 
sibility that the spontaneous fields have almost the same 
magnitude ([ioH e = uje/j^) for all /Lt + sites in the sample. 
Extracted H e (T) data show a similar temperature depen- 
dence as A e (T) and approach 1.9 x 1CT 4 T as T — > (not 
shown) . 

As the origin of the increased relaxation in the SC 
state, the following spurious effects can be ruled out. Any 
inhomogeneous distribution of extremely small residual 
fields (< 1//T < A e (T -> 0)/7 p <C ^ H cl - 5 x 1CT 3 
T caused by the superconducting transition can- 

not explain the observed anomaly. A null experiment 
confirmed that the 4N-silver holder does not exhibit any 
noticeable anomalies. The absence of such an anomaly in 
the previous measurements for UPt3 j^, Celrlns, and 
CeCoIns [i^l, performed under the same conditions, also 
demonstrates that the observed anomaly is uniquely at- 
tributed to the PrOs4Sbi2 specimen. 

Any lattice effects are also ruled out for the following 
reasons: The observed nuclear dipolar field of 1.7 x 10~ 4 
T is close to 2.0 x 10~ 4 T calculated for a most probable 
/i + stopping site (1/2, 0, 0.15) (12e in Wyckoff nota- 
tion, space group Im3), which is one of the interstitial 
sites of Sb-icosahedron cages and is surrounded by six 
nearest Sb ions. This site is the same as the one ex- 
pected for an isostructural PrFe 4 Pi2, as suggested from 
transverse- field ^SR measurements |24j . A simple cal- 
culation shows that a lattice shrinkage of ~ 7% would 
be necessary to explain the observed increase in A when 
A e = 0. This scenario contradicts a thermal expansion 
measurement [25j. Furthermore, the electric field gradi- 
ent at the Sb nuclei 0] does not show any corresponding 
anomalies across the SC transition. 

The spontaneous internal field has two types of possi- 



ble sources depending on the spin and/or orbital parts of 
Cooper pairs having nonzero values: (i) (for both nonzero 
spin and orbital moments) spontaneous undumped su- 
percurrents induced in the vicinity of impurities, sur- 
faces, and/or domain walls between the degenerate SC 
phases, where the order parameter has spatial inhomo- 
geneities |2(| [2t| , and (ii) (for nonzero spin moments) 
a finite hyperfine field induced at the /z + sites. 

An example of the source (i) is found in Sr 2 Ru0 4 [H, 
where A exhibits a spontaneous increase below T c . This 
can be understood by considering that a dilute distri- 
bution of the sources in the sample would result in a 
Lorentzian-type field distribution (2^. Although it is 
possible that the present observation in PrOs4Sbi2 is 
also due to the same type of sources, the spontaneous 
increase appears in A e instead of A. This fact suggests 
that the field distribution is rather of Gaussian type and 
the second moment has a finite value. The mean-free 
path t ~ 2000 A estimated from the Dingle tempera- 
ture (proportional to the electron scattering rate) in the 
dHvA experiments @ , which is much larger than the SC 
coherence length £ = 120 A, suggests that the sample 
is in a clean limit. If the scattering centers determin- 
ing the Dingle temperature are the main sources of the 
spontaneous fields, the average nearest-neighbor distance 
(~ t) would be comparable to the magnetic penetration 
depth A = 3400 A Q. This could explain why the field 
distribution is of Gaussian type. Although quantitative 
comparisons with the models |2y, |22| are hampered by a 
lack of detailed information on the sources, it is certain 
that the field at the center of the sources should be much 
larger than A e /7 M . 

The field source (ii) is possible only if the SC pairing 
symmetry is not only of spin triplet (odd parity), but 
also of "nonunitary" type (analogous to the Ai phase 
of superfluid 3 He 29] ). In such a state, the three- 
dimensional vector function e£(fc), representing the wave- 
vector dependent SC order parameter is complex 
and a nonzero spin moment given by d(k)xd(k)* could 
produce a finite hyperfine field on Note that no 

experimental identification of such nonunitary supercon- 
ductivity has ever been reported. If this is the case in 
PrOs4Sbi2, Mo-ffe corresponds to the "intrinsic" spon- 
taneous field at /i + sites, which is irrelevant to the ex- 
istence of impurities. The magnitude of the field can 
be theoretically estimated to be ~ 7 x 10~ 4 ?y T pjof. 
where rj denotes the degree of the electron-hole asym- 
metry of the conduction band at the Fermi level and is 
of the order of T c /Tp (Tp is the Fermi temperature or 
the Kondo coherence temperature of the heavy quasipar- 
ticle band). A simple estimation of ~ 1 x 10~ 4 T using 
Tp ~ 10 K, which is given by the temperature depen- 
dence of Ti 0]> 1S m good agreement with the present 
observation. A unique feature of the nonunitary state is 
that the spin-up and spin-down Cooper pairs have differ- 
ent sizes of the excitation gap |l5j, |30| . The smaller gap 
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FIG. 3: Longitudinal field Hlf dependence of the relaxation 
rate A. The curves are fits to the model given by Eq. (3). 

branch with <i(A;)| 2 — \d(k)xd(k)*\ ^= might explain 
apparent residual contributions appearing at T <C T c \ in 
some physical quantities [ill fl2| . 

We performed measurements in weak longitudinal 
fields (LF) (hqH^f), parallel to the initial fi + polariza- 
tion, and confirmed that the spontaneous internal fields 
are static on the microsecond time scale. In the mixed 
state after field cooling in /io-ffLF = 100 x 10~ 4 T, the re- 
laxation due to the spontaneous fields is no longer visible 
(already decoupled) , and excellent fits to the spectra are 
obtained using only the exp(— At) term. The tempera- 
ture dependence of A in this field (see Fig- EJ exhibits no 
distinct anomaly through the SC transition. 

In zero field, A shows a gradual increase with decreas- 
ing temperature even in the normal state, in contrast to 
other HF superconductors 0, Efl E3 | . The Hlf depen- 
dence of A shown in Fig. [3] can be nicely fitted by the 
Redfield model: 

A = 2( 7mMo I/ 1oc ) 2 t c /[1 + {^ H lf t c )% (3) 

indicating that this relaxation process is due to fluctuat- 
ing local fields characterized by the magnitude hqHi oc = 
3.7 (4.7) x 10~ 4 T and the correlation time r c = 
0.51 (0.27) /is for 2.0 (4.2) K. We expect that the re- 
laxation process represented by A may reflect the 4/- 
electron dynamics associated with the gap feature given 
by A^cef- 

In conclusion, the ZF-/1SR measurements in PrOs4Sbi2 
have revealed an appreciable increase in the internal mag- 
netic fields below around the onset of superconductivity 
(T c i = 1.82 K). This provides clear evidence for the bro- 
ken TRS in the SC state, which will help us to narrow 
down the number of possibilities for the symmetry of the 
SC order parameter. From the present results, however, 
a definite conclusion cannot be drawn regarding which 
part of the Cooper pairs, spin or orbital (or both), carries 
a nonzero magnetic moment. The broken TRS indicates 



that the SC state belongs to a degenerate representation, 
which has internal degrees of freedom. This is in line with 
the possible existence of multiple superconducting phases 
suggested by the specific heat and thermal conductivity 
studies. Considering the normal state properties, where 
quadrupolar degrees of freedom of Pr ions play a major 
role, the present results point to an exotic conjecture that 
the magnetic and heavy Cooper pairs in PrOs4Sbi2 are 
formed by nonmagnetic interactions mediated by quan- 
tum quadrupole fluctuations. Although this needs to be 
substantiated further by other means, PrOs4Sbi2 pro- 
vides us with a new, and apparently unique, example of 
superconductivity in the field of strongly correlated elec- 
tron systems. 
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